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ABSTRACT 

Turbulent motions induce Doppler shifts of observable emission and absorption lines motivating 
studies of turbulence using precision spectroscopy. We provide the numerical testing of the two most 
promising techniques, Velocity Channel Analysis (VCA) and Velocity Coordinate Spectrum (VCS). 
We obtain an expression for the shot noise that the discretization of the numerical data entails and 
successfully test it. We show that the numerical resolution required for recovering the underlying 
turbulent spectrum from observations depend on the spectral index of velocity fluctuations, which 
makes low resolution testing misleading. We demonstrate numerically that, dealing with absorption 
lines, sampling of turbulence along just a dozen directions provides a high quality spectrum with the 
VCS technique. 
Subject headings: 



Oh; 

6 " 



> 
oo 

o 



oo 

o 



1. INTRODUCTION 

As a rule, astrophysical fluids are turbulent and the 
turbulence is magnetized. This ubiquitous turbulence 
determines the transport properties of the interstellar 
medium (see Elmegreen & Falgarone, [13, Stutzki, [H, 
Balesteros-Peredes et al. [li) and the intracluster medium 
(se e Sunyaev, Norman & Bryan, '39?, Ensslin & Vogt, 
\\%. Lazarian, [20( ). many properties of Solar and stellar 
winds (see Hartman & McGregor, [iGh . and so on. To un- 
derstand heat conduction and the propagation of cosmic 
rays and electromagnetic radiation in different astrophys- 
ical environments it is absolutely essential to understand 
the properties of the underlying magnetized turbulence. 
The fascinating processes of star formation (see McKee & 
Tan, HH, Elmegreen, [91, Mac Low & Klessen, HO) and in- 
terstellar chemistry (see (flSh and references therein) are 
also intimately related to properties of magnetized com- 
pressible turbulence (see reviews by Elmegreen & Scalo, 

We should stress that while the density fluctuations are 
an indirect way of testing turbulence, the most valuable 
information is given by the velocity fleld statistics en- 
coded in spcctrometric observations. The problem with 
those, that has been realized from the very start of the 
research in the field, is that the most^ that one can get 
from such observations are the Position-Position- Velocity 
(PPV) data cubes, where at the every point the image of 
the emitting turbulent volume the Doppler-broaderned 
spectrum is measured. To recover the statistics of tur- 
bulent velocity from the PPV, one has to account for 
the mapping of the turbulent velocity field from the real 
space to PPV data cubes. 

At the moment, several options for turbulence studies 
are available. For instance, using velocity centroids al- 
lows recovering of velocity spectrum (Esquivel & Lazar- 
ian, [13, Esquivel et al. 2007). As shown there, this 

^ In some situations the data is more sparse than that. For 
instance, absorption Unes usually do not provide good spatial sam- 
pUng of the turbulent volume. Moreover, interferometric measure- 
ments of the emitting volumes may not have the entire u-v coverage 
to restore the spatial distribution of PPV emissivity either. 



technique fails to recover the underlying velocity spectra 
for supersonic turbulence. For such turbulence one can 
use the Velocity Channel Analysis (VCA) (see Lazarian 
& Pogosyan, ^o), which deals with channel maps or 2D 
slices of PPV, and Velocity Coordinate Spectrum (VCS) 
(se e Lazarian & Pogosyan, |2l, Lazarian & Pogosyan, 
[ll, Lazarian & Pogosyan 2008, henceforth, LPOO, LP04 
LP06, LP08, respectively, Chepurnov & Lazarian, |3, 
henceforth CL08) which studies the fluctuations of PPV 
intensities along the v-coordinate (see Lazarian 2008 for 
a review). 

VCA has been already used with both atomic hydro- 
gen and molecular line data to recover the underlying 
turbulence spectra (see review by Lazarian, 21). The 
technique has been tested numerically for the optically 
thin case in ((25l ) and (flSl ) and also with absorption effects 
in ((35i ). However, rather strong shot-noise was observed 
in the synthetic channel maps. For instance, numerical 
simulations for VCA in (Ts") show that discretization over 
line of sight may substantially distort the spectrum at 
high wavenumbers. It was also shown for the thin-slice 
(velocity-dominated) case that the impact if shot-noise 
can be minimized by using thickest possible velocity slice. 
The shot-noise was also observed to decrease with the in- 
crease of the cube size used for analysis. Although the 
effect was associated with the discretization of the nu- 
merical data, a quantitative description of the noise was 
missing. Moreover, we feel that this effect confused some 
researchers who started to wonder about the practical 
utility of the technique (32). 

While VCA is an established tool with numerous ex- 
amples of practical applications, VCS is a more recent 
development. Although it was formulated in the same 
paper as VCA, i.e. LPOO, the actual application of VCS 
is just beginning (4). 

In this paper we provide the numerical testing for both 
the VCA and the VCS techniques. In particular, we 
subject the origin of shot noise to scrutiny. 

In in Sect. H we describe general considerations about 
studied turbulent flelds. In Sect. [3l we provide some 
facts about the VCA and VCS. In Sect. HI we analyze 
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shot-noise. In Sect. [5l we describe the resuhs of numer- 
ical testing. The case of observations in absorption hnes 
is discussed in Sect. ID We provide a discussion in Sect. 

El 

2. UNDERLYING FLUCTUATIONS OF VELOCITY AND 
DENSITY 

Our assumptions about underlying turbulence coincide 
with those in LPOO. In particular, in what follows we 
assume that turbulent velocity and density are homo- 
geneous^ and isotropic and admit a statistical descrip- 
tion in terms of correlation (or structure) functions and 
the corresponding power spectra^. The assumption of 
velocity isotropy and its relevance to the actual mag- 
netized astrophysical fluids was discussed at length in 
LPOO, Lazarian & Pogosyan (2004), LP06. There it was 
argued that the degree of anisotropy that is induced by 
the interstellar magnetic field does not compromise the 
power spectrum recovered from observations. This argu- 
ment was supported by the analysis of synthetic maps 
obtained using 3D MHD simulations of turbulence (Es- 
quivel et al. 2003). 

For a vector field like velocity, the correlation function 
and the correspondent power spectrum take the form of 
a tensor. In addition, the power spectrum is different 
for solenoidal and potential field components. For each 
case the tensor properties can be represented by a factor 
that depends only on angles. Additional scalar factor is 
responsible for energy distribution in wavenumber space. 

Because of the self-similarity of turbulence, we expect 
a power law scaling for the power spectrum in the inertial 
range. This power law can be either "long wave domi- 
nated" (or "steep"), if the respective spectral index is 
greater than 3, or "short wave dominated" (or "shallow") 
otherwise (see LPOO). To keep total energy limited, both 
types of spectrum must have a cutoff: at low wavenum- 
bers for the steep one, and at high wavenumbers for the 
shallow one. The cutoff is a physical requirement that is 
satisfied in all astrophysical systems. 

In this paper we assume the velocity spectrum to be 
steep, as a shallow velocity spectrum is not physically 
motivated (see Cho & Lazarian 2005). For the density 
spectrum we consider both steep and shallow spectra. 
The shallow spectrum for density emerges for high Mach 
numbers (see Padoan et al. 2004, Bereznyak, Lazarian & 
Cho 2005, Kim & Rye 2005, Kowal, Lazarian & Beres- 
nyak 2007). 

The quantities we deal with in spectral line observa- 
tions arc the velocity and the gas emissivity. The latter 
can be proportional to the density of gas p (e.g. emission 
lines of HI) or (e.g. recombination lines in plasma). 
The latter regime modifies the analysis. Namely, for 
steep density spectrum (e.g. when ap > 3), the asymp- 
totic slope of the emissivity spectrum is the same, and 
different, — 2ap — 3, for the shallow density spectrum. 

We are interested in velocity statistics, and the gas 
emissivity is considered only because its fluctuations also 

^ Random fields are homogeneous in the statistical sense when 
their statistical properties do not depend on the position. This al- 
lows us to keep only the dependencies of our correlation functions 
on the distance between the points. Homogeneity is also a require- 
ment for providing the volume averaging that observers would rely 
on to apply our results 

^ The latter pair of functions is related to each other through 
Fourier transform (see Monin & Yaglom 1975). 
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Fig. 1 . — An illustration of the mapping from the real space to the 
PPV space. In the real space 3 eddies above have the same spatial 
size, but different velocities. They are being mapped to the PPV 
space and there they have the same PP dimensions, but a different 
V-size. The larger is the velocity of eddies, the larger the V-extend 
of the eddies, the less density of atoms over the image of the eddy. 
For the eddy 3 the slice is "thin", i.e. the statistics is defined 
by the velocity spectrum. The images of eddies corresponding to 
velocities less than the channel thickness is different. The velocity 
slice is "thick" for the eddy 1, i.e. the statistics is defined by the 
emissivity spectrum, and the eddy 2 is in an intermediate regime. 

affect the observable PPV fluctuations. The influence of 
the fluctuations of gas emissivity is mitigated when they 
have a steep spectrum. In this situation, in fact, they 
could be ignored within the VCS technique. However, 
when fluctuations of gas emissivity follow shallow spec- 
trum these fluctuations should be properly accounted for. 

3. BASICS OF VCA & VCS 

Fig. [1] illustrates the essence of the VCA technique. 
The data in the PPV data cube are analyzed using slices 
of AV^ thickness. This corresponds to the analysis of in- 
tensity fluctuations in channel maps (see Green 1993), 
which predates the technique. The gist of the technique 
is a relation between the variations of the spectral in- 
dex of channel map fluctuations and the thickness of the 
channel maps. The differences in the statistics of the 
fluctuations stem from the fact that the images of the 
eddies can have an extent both larger and smaller than 
PPV slice thickness AF. The former case is termed by 
LPOO "thin slice" , and the latter is termed "thick slice" . 

In Fig. [21 we illustrate the essence of the VCS tech- 
nique. Depending on the resolution of the telescope one 
may or may not resolve the spatial extent of the eddies 
under study. This also results in two distinct regimes of 
turbulence studies. The ability to study turbulence in 
the turbulent volumes that are not spatially resolved is 
a unique ability of VCS. 

Below we briefly overview the main analytical results, 
obtained in LPOO, LP04, LP06 and CL08, which are rel- 
evant for this paper. We, however, do not deal with the 
case of VCA/ VCS studies of turbulent volumes where 
self-absorption is important"*^ (LP04) or the case of sat- 

* In this case analysis is more simple if we use Fourier space. 
Absorption effects can be described only in real spa ce, wh ich im- 
plies the use of structure functions instead (see | |24| ). II23I) ). It is 
shown in |24i) that for VCA absorption can be neglected as soon as 
the scale of turbulence under study is less than the physical length 
corresponding to the optical depth of unity. A similar criterion is 
applicable to VCS dH). 
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Fig. 2. — VCS technique: effects of resolution. The fluctuations 
along the velocity coordinate are analyzed. Eddies within the tele- 
scope size beam, e.g. eddy 1, are in a low resolution mode. Eddies 
with the size exceeding the one of the beam, e.g. eddy 2, are in 
the high resolution mode. 

urated absorption lines (LP08). In terms of the presen- 
tation of our results, below, we use the approach and 
notations from CL08. 

Our calculations are based on the following expression 
for a spectral line signal, measured at velocity vq and 
given beam position e (see Appendix |A] for the details) : 

S{e,vo)= j w{e,v)dve{v)f{vr{v)+v:-a{v)~vo). (1) 

where e is normalized emissivity, v^'^^ is a line-of-sight 
component of the regular velocity (e.g. the velocity aris- 
ing from the galactic velocity shift), Vr is the line-of- 
sight component of the random turbulent velocity, / de- 
notes the convolution between channel sensitivity func- 
tion^ and Maxwellian distribution of velocities of gas par- 
ticles, defined by temperature of emitting medium, and 
w is a window function defined as follows: 



(2) 



where Wb is an instrument beam, which depends on an- 
gular coordinates </> and 9, while We is a window function 
defining the extent of the observed object. For instance, 
for a turbulent volume confined by two planes perpen- 
dicular to the line of sight, the window function is as 
follows: 



We{z) 



1, 
0, 



z e 



_zoiZi\ 
'zo,zi] 



(3) 



In what follows, we consider the Fourier transform of 
a spectral line^: 



= f{kv) J w{e,r) dr- 

£(r) exp(-z/ct,(wr(r) -I- v^^^^r))). 



(4) 



This function can be easily determined from observa- 
tional data. 

I.e. amplitude-frequency response of the channel normalized 
to the integral value with frequency in velocity units 

Variable fc„ plays here the role of in LPOO, being, however, 
different in dimension [kz = bkv, see Eq. l(7)l. We use it here to 
avoid complications when b = 0. 



If we correlate S taken in two directions 1 and 2, we get 
the following measure, which can be used as a starting 
point for the mathematical formulation of both the VGA 
and VCS techniques (see CL08 for details): 



i^(ei,e2,A:„) = (^S{ei,ky)S*{e2,ky) 

= P{ky) J w{ei,r) dr J w{e2,r') dr' ■ 
(e(r)e(r')) {exp{-iK{vr{r) - «.'(r')))> ' 
exp(-iA:„(«7f(r)-<f3(r'))), 

(5) 

The first averaging gives us an emissivity correlation 
function Ce(r — r'). Averaging of the exponent can be 
performed with the assumption that the velocity statis- 
tics is Gaussian^: 



{exp{-iky{vr{r) - Wr'(r')))) 
= exp[-^^{ivrir)-v.Ar')r) 



(6) 



To proceed, we assume that the beam separation and the 
beam width are both small enough that we can neglect 
the difference between Vr and Vz (we consider z-axis to 
be a bisector of the angle between beams). We also as- 
sume that u^*^^ (r) depends only on z and admits a linear 
approximation: 



where b characterizes regular velocity shear. 
This leads us to 



(7) 



/^(ei,e2,fc„) = Pikv)- 

J w{ei,r)dr f w{e2,r') dr' C^{r ~ r')- 

exp (^~!^Dyz{r - r') - ikyb{z ~ z')^ , 

(8) 

Having substituted r — r', we can write the following 
expression for K: 



K{ei,e2,ky) ^ P{kv) / .g(ei, 62, r) dr • 



Ce(r) exp {-^D^zi^) - ik^bz^ 



where 



g(ei,e2,r) = J 'w{ei,r')w{e2,r' + r) dr' , 



(9) 



(10) 



/ is a Fourier transform of an effective channel sensitivity 
function /, Cg is an emissivity correlation function, and 
D-uz is a velocity structure tensor projection: 



D,z{r-r')EE{{vzir)-Vz{r'))') 



(11) 



If we set g — (j(R — 'Ri2)w^^a{z), where R12 is a 
beam separation in pictorial plane^, Wb^a and We^a are 

^ We assume that the velocity field has a Gaussian Probability 
Distribution Function (PDF). The latter corresponds well to both 
experimental (Monin & Yaglom 1976) and numerical (Biskamp 
2003) data. 

* The pictorial plane is the plane perpendicular to the line of 
sight. It provides us the coordinates in spatial units for the beam, if 
we consider beam's cross-section by a pictorial plane at the distance 
of a remote object. 



4 



Chepurnov & Lazarian 



auto-convolutions of Wh and w^, which corresponds to 
the case of remote emitting object, we have for K: 

i^(Rl2, ky) = PiK) J WtM^R - Rl2)We,aiz) dv ■ 

Ce(r) exp {-^D^^{v) - ik^bz^ 

= fHk,)Jwl{K)dK- (12) 
/iye,a(2)dre-'K(R-Ri.). 

Ce(r) exp (^~^Dy^{r) - 

Taking a Fourier transform over R12, we have the power 
spectrum, designated in VGA as P^: 

Ps{K,ky) ^p{k,)wt\K)Jw,^,dv 

g-.(KR+6fc„.)c.^(j.)exp(-f i^„,(r)) , ^^^> 

which coincides with the corresponding expression in 
LPOO. 

It is convenient to introduce P3 as a 3D power spec- 
trum of a PPV data cube with assumption of infinite 
resolution over v and R: 



(14) 



P3(K,A:„) = (|^(K,fc,)|2^ = /u;,,„dr• 
g-.(KR+6fc„^)c^(J.) exp (-^i^„,(r)^ 

In this case it is easy to see that Ps(K, fc^,) — 

/2(fc„)5^'(K)P3(K,fc„). ^ 

VGA studies the two-dimensional spectrum of fluctu- 
ations within a velocity slice of a PPV cube: 

P2(K) ^ ^°^^dk,Ps{K,k,) 

= G^'(K)/_!°^ dfc„/2(fc„)P3(K,fc,), ^ > 

which can be easily determined from observations. This 
technique provides a way to determine the slopes of ve- 
locity and density power spectra, Uy and a^. In fact, for 
a channel with small enough effective width^ (see Fig. 
[1]), P2 is in the velocity dominated regime with slope 
(9 — ai,)/2, provided that the density has a steep spec- 
trum, i.e. tte > 3. 

Whether or not the latter is true can be established 
through using from column density maps. As we neglect 
the effects of self-absorption, the column densities can be 
obtained via v-integration of PPV cubes. Naturally, in 
the column density maps the spectrum is affected only 
by density and its slope is a^. The situation is a bit 
more complicated when the density is shallow (i.e. when 
< 3) , which is a case of high Mach number turbulence 
(see Beresnyak, Lazarian & Gho 2006) and the density 
combines with velocity to affect the fluctuations in thin 
channels. For a more detailed analysis, see LPOO. The 
predictions for different cases in the velocity-dominated 
regime are summarized in Table [T] Note that these re- 
sults are obtained by assuming that the turbulent object 
under study is so distant that the geometry of parallel 
lines of sight is applicable. 

The one-dimensional spectrum Pi is the subject of the 
VCS studies. It corresponds to the case when the two 
beams involved in K coincide, i.e. 



Pi(fc,) =i^(Ri2,fc,)|R,,=o = /'(fc.) / g{r)dv 
Ce{r) exp (^-!^Dyz{r) - ikybz\ 



(16) 



TABLE 1 

VGA PREDICTIONS ABOUT P2 SPECTRAL INDEX, 
STEEP DENSITY 



density spectrum 


2-d spectrum 


1-d spectrum 


steep 


9— ci„ 
2 


7-c«„ 
2 


shallow 


2ae-a„+3 






2 


2 



Note. — Here av is the velocity spectral index, 
is the density spectral index. Density spectrum 
is considered to be "steep" if > 3 and "shallow" 
otherwise. 



TABLE 2 

VCS PREDICTIONS ABOUT Pi SPECTRAL INDEX, PARALLEL LINES 
OF SIGHT 



density spectrum 


pencil beam 


flat beam 


low resolution 


steep 
shallow 


2 


4 


6 


a^—Z 
2(01^-2) 


a„-3 
2(c^-l) 


a„-3 
2ap 


a„-3 


a„-3 





Note. — Here is the velocity spectral index, is the 
density spectral index. Density spectrum is considered to be 
"steep" if > 3 and "shallow" otherwise. 

which can be written also in terms of P3 as: 



Pi(fc,) = /"(fc.) / dK5);^(K)P3(K,fc,). 



(17) 



® It is approximately equal to y/iSv^^ + {2vt)^), where Sv^^ is 
the channel width and vt is the thermal velocity of emitting atoms. 



As it was shown in (f23h and (0), here we also have two 
spectral regimes. They depend on beamwidth (see Fig. [5] 
for illustration). For the high resolution mode {k^ is less 
than the velocity variance on the beam scale) the slope 
of Pi is 2/{av — 3), otherwise, in the low resolution mode 
it is 6/(q!^ — 3) (a steep density spectrum is assumed for 
the both cases). A more complete list of the predictions 
for Pi is presented in Table [2l The flat beam there cor- 
responds to the telescope beam for which the resolution 
along one of the spatial axes is infinitely better than the 
resolution for the perpendicular axis. 

If we compare Eqs. (fT5|) and (fT7|) . we see that the VGA 
and the VGS provide complimentary approaches and are 
related to the same quantity, P3. Their properties must 
therefore be closely related. 

4. SHOT-NOISE ARISING FROM NUMERICAL FINITE 
RESOLUTION 

Earlier attempts to test the accuracy of VGA faced a 
problem; namely the distortion of the spectrum of fluctu- 
ations within a v-slice of the PPV data cube (see Lazar- 
ian et al. 2001). The effect from the very beginning was 
associated with the discretization of the numerical data, 
but attempts to remove the noise by the interpolation of 
the numerical data were not successful. In particular, in 
(jUt ) it was shown that the quality of the spectrum in the 
v-slice depends on the number of points along the line 
of sight in the numerical cube used to create the PPV 
cube. The interference noise was termed shot-noise to 
reflect its stochastic nature: insufficient statistics for the 
signal in spectrometer channels^*'. Below we show that 
the same type of noise is present in the VGS analysis. In 

If we take too low Nz in our simulations, there can be too 
few points that contribute to the signal in a particular channel, 



which makes such a value "unreliable" 



However, one can face 
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this case, Pi asymptotically approaches a constant min- 
imum value for large fct,, which diminishes the spatial 
extent over which one can study turbulence (see Fig. [71 
right). 

To have a realistic spectral line we need a velocity field 
to be adequately resolved over line of sight (see also In- 
ogamov & Sunyaev, |T3). This means that the typical 
scale over z for which velocity stays within a single chan- 
nel near its extremum, should have enough discretization 
points, because such extrema give the dominant contri- 
bution to a spectrometer signal. 

Let us discuss the required resolution, related to the 
size Nz of the data cube along z-axis. One can con- 
sider the non-linear space- velocity mapping^^ induced by 
the turbulent velocity. The relation between the scale of 
real space I and the corresponding velocity scale vi is 
given by vi ~ ;("f-3)/2^ which for the Kolmogorov tur- 
bulence, i.e. a — 11/3, gives the familiar vi ~ /^/^ rela- 
tion. The minimal scale Imin which can be studied with 
either VCS or VGA technique corresponds to the dis- 
tance between adjacent points in the original data cube. 
This minimal distance translates into the smallest veloc- 
ity over which the velocity fluctuations can be studied, 

I.e. Vz.min 



(a.-3)/2 



This minimal velocity corre- 



sponds to the maximal wavenumber along the velocity 
coordinate, which can be recovered from the simulations 

kv,m,ax ^ ^/^z,rnin- Therefore ^v^maa: ^ ^min ' the 

data cube has Nz points, Imin ^ ^7^- This results in 

kv,max ~ Nz" For kv numbers to contain the infor- 

mation about turbulence ky should be less than ky^„iax- 
Consequently the number of points along the line of sight 

in the numerical data cube should scale as kv^max 
This provides the same dependence as the more detailed 
Eq. dill). 

Naturally no VGA studies are meaningful below the 
channel thickness of l/fcu.maa;- 

The arguments above also mean that, according to Ta- 
ble 2, the required number of points for flat beam scales 

as kt(max and for low resolution as ky^mcu: ^\ which 
provides even steeper dependencies than the pencil beam. 
Thus, testing our formulae with low resolution will be 
more challenging. 

To estimate the exact resolution dependence, we have 
to relate the correspondent typical scale over line of sight 
with the channel width. The natural choice is to calcu- 
late the velocity dispersion for wavelengths shorter than 
this scale. As mentioned above, the contribution of the 
velocity profile shapes other than extrema, provided by 
the remaining lower harmonics, can be neglected. There- 
fore, the scale found by reverting this expression is what 
we seek. 

In other words, we divide the ID velocity spectrum 
domain into the intervals separated by the scale in ques- 
tion La, and, as spoken, assume that the harmonics in 
its low- wavenumber part provide an extremum of veloc- 
ity. For the higher harmonics, in order to place the bulk 
of the signal inside one channel, we should demand that 

similar situation in X-ray observations, where individual photons 
are counted. In this case the time of observation plays the role of 

Such mapping can be illustrated by an example provided on 

Fig.m 



the correspondent variance g^La) is about half a channel 
width: 



<j{La) 



(t(oo) 
2iV„ ' 



(18) 



where La is a scale at which velocity stays within a single 
channel near its extremum, N^, is the number of chan- 
nels per velocity variance, and <j{L) is a velocity variance 
for scales smaller than L, expressed through ID power 
spectrum^^ E{k): 



a^{L) = 2 / E{k)dk 



(19) 



If a 3D power spectrum has a cutoff at fco ~ 27r/i, 
where L is an injection scale, E{k) has the following form 
for the isotropic 3D spectrum of Vz 



E{k) 



1 



Uv — 2 



/c-(""-2), |fc| > fco 



\k\ < fcn 



Then, having solved Eq. for La we get 



La 



2N„ 



(20) 



(21) 



If UA is a number of points per La, then the total 
number of points along z is as follows: 



Nz 



UA- 



2N„ 



\Jay-2 



(22) 



where L^ is the size of emitting structure, ua is deter- 
mined from simulations and is experimentally found to 
be around 6. 

The requirements for Nz are rather demanding. If we 
take our typical settings: 128 channels per triple line 
width and LsjL ~ 2, then for — 4, Nz should be 
at least 1.1 • 10^ points, for = 11/3 the corresponding 
number is 4.3-10^, while for a„ = 3.5 the number is as big 
as 1.7 • 10^. Therefore it is not surprising that without 
this knowledge some researchers could get puzzled by 
their results (e.g. Miville-Deschnes et al. |32| ). 

Gondition (I22|) guarantees that the spectrum is clear 
from shot-noise for all available fc^'s, up to 27r/Ai;, where 
Aw is the channel width. This estimation refers to the 
case without smoothing along angular coordinates, i.e. 
for the high-resolution regime. The low-resolution VGS 
spectrum drops faster and therefore the requirements are 
even more restrictive. 

5. RESULTS OF NUMERICAL TESTING 

5.1. Recovering Spectral Indexes 

VGS spectra, illustrating the condition given by 
Eq. ([22]) are presented in Fig. [5l The velocity spectral 
index decreases from top to bottom (the values are 4, 
3.78 and 3.56). According to our estimates above, this re- 
quires a drastic increase of the number of points. Indeed, 
Nz is, respectively, 2^**, 2^^ and 2^^. In the left column, 
with shot-noise clearly present, Nz is 4 times underval- 
ued. Examples of restoration of the velocity spectral in- 
dices from simulations are presented in Tab. [3l They 

This Should not be confused with Pi, which refers to PPV 

space 
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TABLE 3 

Restoring of velocity spectral 

INDEX av from simulations 



expected a„ 


restored ol^ 


Nz 


3.56 


3.54±0.01 


8388608 


3.67 


3.65±0.02 


524288 


3.78 


3.79±0.02 


131072 


3.89 


3.89±0.02 


32768 


4.00 


3.91±0.04 


16384 




Cluster 
of stars 



Note. — 16 realizations of spectral 
line have been used for each Pi 

show good agreement. Values of have been chosen 
regarding Eq. (j22p . Thus, knowing the requirements on 
A^2, we can generate the VCS and VGA spectra (see Fig. 
H with a„ = 4). 

But for shallower velocity spectra, satisfying such con- 
ditions for VGA may be difficult, even if we use 2-d data 
instead of a 3-d data cube. For example, the condition 
for the VGA simulations at a.^ = 11/3 implies a 
square array with side 420000, which is difficult to im- 
plement at present. If we decrease this size, the VGA 
spectrum gets distorted as seen in Fig. [71 left. Therefore 
instead of prescribing the number of VGA channels we 
have to seek alternative approaches. 

5.2. Choosing Optimal Thickness of Slices 

In (flsl ) it was shown, that if we "optimize" the thick- 
ness of the velocity slice (which is equivalent to chang- 
ing the velocity resolution), we get a correct spectrum. 
Let us consider why this can happen. Looking at the 
right panel on Fig. [71 we see that at higher ky the VGS 
spectrum is affected by shot-noise. As both VGA and 
VGS spectra are just different projections of the same 3- 
d power spectrum in the PPV space (see Eqs. (HSj) and 
(fT7)) ). the noisy harmonics seen in the Pi plot are most 
likely to affect the P2 spectrum. If we use only fc„'s, 
which correspond to the intact part of Pi, the picture 
improves significantly, as seen in Fig. [51 left. Such trun- 
cation of the fc.u-domain is equivalent to taking a thicker 
velocity slice in (fl3| ). Rather thick slices were used in a 
numerical study by Padoan et al. (2006), which explains 
why they did not face the problems related to shot-noise. 

It is also important to test the VGS prediction of steep- 
ening of Pi when the beam is made wide enough to emu- 
late the low resolution regime. To do so we used steeper 
velocity spectrum, with slope 4, which is less demanding 
in terms of the number of points along line of sight. If we 
meet condition (|^^ we have a good agreement with the 
predictions for both the high- and low-resolution regimes 
(see Fig. [ID. 

6. STUDIES OF TURBULENCE WITH ABSORPTION LINES 
FROM POINT SOURCES 

As clear from the earlier work (see LPOO), as well as 
from our presentation above, the VGA and VGS are re- 
lated techniques. We must stress that although fields of 
applicability of discussed techniques intersect, they do 
not coincide. For instance, VGA requires less resolution 
over the velocity coordinate, because the range of scales 
can be obtained with spatial variations. Thus, poten- 
tially, it can be used for studies of gas with Mach number 
lower than what is required for VGS. 



Absorbing 
medium 




Effective tjeam 



Fig. 3. — VCS allows recovery of velocity statistics from ab- 
sorption lines from stars. The whole Pi is guaranteed to be in 
high-resolution mode, which provides better dynamical range over 
ky. Numerical simulations show, that very few independent mea- 
surements are needed to gain required statistics (see Fig. I10|l . In 
the case of absorption against an extended source we have situa- 
tion like shown on Fig. [2]for an emission line: the eddies can be in 
low (1), high (2) or intermediate (3) resolution mode. The effec- 
tive beam is defined by the background object shape in this case, 
if it isn't resolved by the telescope. For absorption lines observed 
against CMB, we have the picture identical to shown on Fig. [2] 

On the other hand, even if we have very poor statistics 
over angular coordinates and VGA does not't work, VGS 
still can be used. This opens new possibilities for studies 
of turbulence with absorption lines. In general, for ab- 
sorption with extended sources both techniques can be 
readily applicable. As the spatial scale of the source de- 
creases, the utility of the VGA decreases as well. Indeed, 
a lot of point-like sources are required to provide the 
spatial statistics required for the VGA. On the contrary, 
VGS utilizes sampling over different directions just to get 
different statistical realizations of the turbulent process. 
Thus the number of point-like sources (e.g. stars) may 
be substantially reduced. The typical observational sit- 
uation for probing turbulence with emission from point 
sources is illustrated in Fig. [3] What is the minimal num- 
ber of sources of absorbed emission that are required for 
restoring the spectrum of turbulence between the sources 
and the observer? Below we answer this question using 
numerical simulations. 

A general case of VGA studies with absorption lines is 
discussed in Lazarian & Pogosyan (2008). There the cri- 
terium for restoring the turbulence signal from saturated 
absorption lines was obtained. For the sake of simplic- 
ity, here we discuss unsaturated absorption lines, arising 
from the medium with total optical depth less than unity. 
While, in general, a number of spectral lines from a sin- 
gle source may be used, we limit our study to the situ- 
ation of a single absorption line. Our simulations show 
that only 10 independent measurements are needed to 
gain required statistics for « 20 (see Fig. [TOj). This 
opens avenues for studies of turbulence with background 
star even when the number of stars behind the turbulent 
volume is limited. No other presently known technique 
other than VGS can recover the velocity spectrum in this 
situation. 

We may note that the number of absorption sources 
may be even fewer than we mentioned above. This is the 
case, for instance, when the turbulent volume sampled 
along the line of sight contains subvolumes with differ- 
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Improving statistics by regular velocity sliear 




One real cation, specmjm compresston is 16. 
Avdue to galaxy rotation is 8 times greater 
tnan turbulent line wkdtn. 



Fig. 4. — VCS can be applied even to a single spectral line. 
Configuration can be as shown on this picture. Regular velocity 
shift from the galactic rotation results in statistics being sufficient 
to recover. Pi 

ent regular velocities. Such a velocity shift is possible 
due to galactic rotation, for instance. In this situation 
the limiting case for study is studying turbulence sam- 
pled by an absorption line from a point source along a 
single direction. In the presence of regular velocity differ- 
ence between different turbulent volumes, the statistics 
may be rich enough to enable averaging of adjacent ky- 
harmonics (see Fig. 2]). 

Naturally, when one uses point-like sources to test tur- 
bulence with absorption lines, the effective resolution for 
the VGA is infinite. In the case of absorption lines aris- 
ing from sources of finite angular size, the situation may 
be different. In particular, we may be also in the low 
resolution regime (see Fig [3]). 

7. DISCUSSION 

Our study in the paper above clarifies a number of is- 
sues related to the studies of turbulence with the VGA 
and VCS techniques. In particular, our numerical testing 
proves the analytical expressions which are corner stones 
of the VGA and VGS techniques. We also identified the 
source of shot noise reported in the earlier numerical test- 
ing of the VGA, provided an analytical estimate of the 



velocity slice thickness required to eliminate the noise, 
and checked successfully that when the requirement of 
the velocity thickness is satisfied, the velocity spectrum 
can be successfully recovered. 

In our treatment we stressed that the VGA and VGS 
techniques are closely related (cf. Eqs. (fT5|) and p7|) ) 
and therefore the presence of shot-noise in VGS in- 
evitably affects the quality of VGA spectrum. We showed 
that the requirement given by Eq. (|22p to eliminate the 
shot-noise in synthetic PPV cubes is rather demanding as 
the velocity spectrum gets less steep, which send a warn- 
ing signal to brute-force attempts to test the techniques 
for an arbitrary index of turbulence spectrum. We con- 
firm, however, an earlier claim in Esquivel et al. (2003) 
that the shot-noise interference is, exclusively, a prob- 
lem of synthetic observations, rather than the real ones. 
Shot-noise will appear in real observations too, when the 
number of received photons is low enough. Incidentally, 
the condition (|22p gives in this case the required number 
of photons and can be used for the estimation of obser- 
vational time required. 

The limitation of our present study was that we dealt 
only with testing of the asymptotical solutions for the 
VGA and VGS techniques. Real data with its limited 
dynamical range of observed velocity fluctuations may 
benefit from fitting the observations of the integral ex- 
pressions rather than the asymptotical solutions. Indeed, 
the actual spectrum in the velocity slice may be thin for 
large eddies and thick for small eddies. Therefore asymp- 
totical solutions may fit only part of the actual spectrum 
of the slice, which will deviate from a power-law. We 
shall deal with the numerical testings of these situations 
elsewhere. 
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APPENDIX 
SPECTRAL LINE 
Let us write out the space-velocity distribution of density: 
no(r, v) = no{r)(p{v - Vmacroir)), 
where (p{v) is Maxwellian distribution of the emitting atoms (ions): 
1 



: exp 



ks is Boltzmann's constant and rio(r) is an integral density. 
In these terms emissivity can be written as follows: 

£o(r, v) = eo{v)ip{v - Wmacro(l")), 

where 

eo(r) = 7ino(r) 

- for linear emissivity law, or 

eo(r) = l2nl{Y) 

— for quadratic emissivity law. 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 
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Here 71, 72 are the coefficients depending on the correspondent cross-sections, and v defines emission frequency shift. 
Let us write out the elementary power^'^ , radiated by the elementary volume dr dv into the solid angle dfl : 



(A6) 



dP ^eoir,v)^dr dv 

= ji:^£oi'r)ip{v - Vfiow) ds dr dv 

where ds is an elementary area at the distance r from the source. 
Then we can express the elementary signal at the spectrometer output as follows: 

1 A tt dP r 

dS=:r, "-^ -Mv-vo)wbA0,^)drdv (A7) 

2fc_B ds /o 

where Aef / is effective area of the instrument, fs (v) is the amplitude- frequency response of the spectrometer normalized 
to the integral value, /o and vq are central frequency and velocity of the channel and Wb,i{0, (j)) is the instrument beam 
normalized to the maximal value. S is measured in Kelvins. 
Accounting for 

Aeff^Sa^ (A8) 

^a-- J wbAe,^)dn (A9) 

where fla is the beam solid angle, Aq is central wavelength and Sa is the aperture efficiency and defining 

" ^ (AlO) 
we obtain: 

1 

dS ^- — -(p{v-Vfiow{r))fs{v~VQ)wb{e,r)dr dv (All) 
Zkb 47rr^ 

where e is the beam direction. Integrating dS we finally have the emission line profile: 

^(e, vo) = f dr !fi^i%^eo(r)/(t;o - Vfi^^r)) (A12) 

where effective channel sensitivity function / is defined as follows: 

00 

f{v)= [ ip{v' + v)f,{v')dv' (A13) 



Random field £o(r) is not homogeneous, at least because the emitting structure is limited in space. To model this 
we introduce a homogeneous field e(r) and a deterministic factor We{r) setting up the borders of an observed object. 
We will also "pack" into e(r) all constant factors: 

w,{r)s{r) ^ ^eo(r). (A14) 
Introducing a window function w as follows^"*: 

w(e, r) EE -^Wfc(e,f)we(r), (A15) 
we finally obtain: 

5(e,«o)= j w{e,r)dre{r)f{vr{r)+vl-a{r)~vo). (AI6) 

CODE 

The simulation programs arc written in CH — h in single-processor and MPI variants. The source of simulation 
programs is available at www. astro. wise. edu/^lazaricin/simulations, see Tab. |4]for the complete list. 



we neglect the self-absorption effects 
^"^ This composite window function is a product of the 3D instrument sensitivity function and the window function defining the object 
extent. 
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TABLE 4 
List of simulation programs 



program 


dimensions 


techniques 


architecture 


spectld 


1 


vcs 


single processor 


spect2d 


2 


vcs 


single processor 


spectSd 


3 


vcs 


single processor 


p_spect2d 


2 


vcs 


MPI 


p_spect3d 


3 


vcs 


MPI 


vca3d 


3 


VCA 


single processor 


p_vca2d 


2 


VGA 


MPI 


p_vca3d 


3 


VCA 


MPI 



Note. — The programs are available at 
www . astro .wise. edu/~lazarian/simulat ions 
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Fig. 5. — Illustration to the requirement |22j to the number of points over line of sight. The rows correspond to different velocity spectral 
index Ot, (4, 3.78 and 3.56 from top to bottom). VCS spectra in the middle column correspond to Nz close to the estimation (2^^, 2^"^ and 
2^^ respectively), in the left column Nz is 4 times lower, in the right — 2 times greater. The cases with low A^^ exhibit shot-noise, which 
results in points getting above the solid line, which represents the analytical expectation. 




1.0 1.5 2.0 2.5 5.0 -5.5 -5.0 -2.5 -2.0 -1.5 -1.0 

og K log kv 



Fig. 6. — Simulations with the number of points over line of sight Nz = 32768, sufficient to suppress shot-noise {Nz = 20000 is required). 
Left: VGA spatial spectrum (2D simulation), shallower solid line - expected velocity-dominated spectrum, steeper solid line — density- 
dominated spectrum. Right: correspondent VCS spectrum for high resolution mode (ID simulation), solid line - expected slope. Velocity 
spectral index is 4. 
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1.0 l.b 2,0 2.i 0,0 -0,0 -2,5 -2,0 -1,5 -1,0 

K log kv 



Fig. 7. — Simulations with insufficient number of points over the line of sight (N^ = 32768 instead of ~ 420000). Left: VGA spatial 
spectrum (2D simulation), shallower solid line - expected velocity-dominated spectrum, steeper solid line - expected density-dominated 
spectrum; the spectrum is distorted (the velocity-dominated one is expected). Right: correspondent VCS spectrum for high resolution 
mode (ID simulation), solid line — expected slope, shot-noise clearly visible. Velocity spectral index is 11/3, which implies larger minimal 
Nz, than for the case on Fig. \E\ 



P2 for dv=375 




Fig. 8. — The same as Fig. but with noisy fci,-harmonics filtered out when calculating the VGA spectrum (left). 




Fig. 9. — 2D VCS simulations for high-resolution (left) and low-resolution (right) modes. Solid lines show expected slopes. Velocity 
spectral index is 4. 
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Fig. 10. — VCS spectrum for 10 realizations in high-resolution mode. This illustrates the opportunity to recover velocity spectrum from 
absorption lines against a limited number of stars. 



